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Abstract 

Analytic formulae are given for QCD corrections to the lepton spectra in decays of 
polarised up and down type heavy quarks. These formulae are much simpler than 
the published ones for the corrections to the spectra of charged leptons originating 
from the decays of unpolarised quarks and polarised up type quarks. Distributions 
of leptons in semileptonic A c and A& decays can be used as spin analysers for the 
corresponding heavy quarks. Thus our results can be applied to the decays of 
polarised charm and bottom quarks. For the charged leptons the corrections to 
the asymmetries are found to be small in charm decays whereas for bottom decays 
they exhibit a non-trivial dependence on the energy of the charged lepton. Short 
life-time enables polarisation studies for the top quark. Our results are directly 
applicable for processes involving polarised top quarks. 
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1 Introduction 



Inclusive semileptonic decays of heavy flavours play important role in present 
day particle physics. In near future with increasing statistics at LEP and 
good prospects for £?-factories quantitative description of these processes may 
offer the most interesting tests of the standard quantum theory of particles. 
At the high energy frontier semileptonic decays of the top quark will be 
instrumental in establishing its properties Jl|. |2|] . 

Recent theoretical developments, discovery of Heavy Quark Symmetries 
[§, H and Effective Theory || as well as subsequent applications of an expan- 
sion in inverse powers of the heavy quark mass uiq || , have led to a consistent 
treatment of semileptonic decays of charmed and beautiful hadrons. In the 
limit vtlq —>■ oo the inclusive decays of heavy flavour hadrons are described 
by the decays of the corresponding heavy quarks. For finite quark masses 
there arise non-perturbative corrections due to the heavy quark motion in- 
side the heavy hadrons |7|]. A systematic approach to these corrections based 
on l/rriQ expansion is described in || and other articles cited in ||. In the 
present article we consider only perturbative QCD corrections to the decays 
of heavy quarks. Non-perturbative effects have to be treated according to the 
recipies provided by the simple model of 0j or by the more formal approach 
of || . On the other hand inclusion of perturbative corrections is necessary for 
a quantitative description of semileptonic decays because these corrections 
are of the order of 20 percent. 

Interesting new opportunities are provided by decays of polarised heavy 
quarks. It is well known [p], that the heavy quarks produced in Z° de- 
cays are polarised. According to the Standard Model the degree of longi- 
tudinal polarisation is fairly large, amounting to (Pb) = —0.94 for b and 
(P c ) = —0.68 for c quarks 0. The polarisations depend weakly on the pro- 
duction angle. QCD corrections to the production cross section have been 
calculated recently ||. Unfortunately, this original high polarisation is to 
large extend washed out during hadronisation into mesons. At the time be- 
ing only charmed and beautiful A baryons seem to offer a practical method 
to measure the polarisation of the corresponding heavy quarks. This method 
proposed by Bjorken || has been successfully applied to both semileptonic 
|i~0|l and hadronic Jill decays of Lambda baryons. The polarisation transfer 
from a heavy quark Q to the corresponding Aq baryon is 100% |12[ , at least 
in the limit wiq — > oo. In view of the growing sample of A c and A^ baryons 
produced at LEP there arises an opportunity to measure the polarisation of 
the c and b quarks originating from the decays of the Z bosons. The angular 
distributions of charged leptons [|13|, [14| and neutrinos fl5j from semilep- 



tonic decays of A c and A b can be used as spin analysers for the decaying 
heavy quarks. Polarisation studies are much more direct for the top quark 
which is a short lived object and decays before hadronisation takes place. 
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Angular distributions of leptons from semileptonic decays of the top quark 
can provide a very interesting information on the space-time structure of the 
corresponding electroweak current 



TBI 



First order QCD corrections to the energy spectra of charged leptons in 



the decays of c and b quarks have been calculated analytically in ref. HT7 
and |Tj|. The results of a later analytical calculation presented in ref. [[19 
disagree with those of jl7| and 
and bottom 



and agree with a Monte Carlo for charm 
20j as well as with numerical results for top |2T|. In |H] 



the 

formulae from |2(J were used and the accuracy achieved was good enough 
to observe discrepancies with [O] and IJ|. The formulae for the virtual 
corrections used in p0[ and |ZT| were taken from the classic articles on muon 
decays [E3]. 



The formulae given in the present article agree with the ones given in [19 



but are much simpler. They agree also with those given in |23j for the joint 
angular and energy charged lepton distribution in top quark decays. In the 
latter case simplification is even more striking. Some numerical results for 
polarised charm and bottom quarks have been already presented in [15]. 

In section |2| the formulae are given for the QCD corrected distributions of 
leptons originating from weak decays of polarised heavy quarks: subsection 
|2.1| contains our main result - the triple differential angular and energy dis- 
tributions of leptons, and in subsections [2.2| and [2.3| the double differential 



distributions are described. The corresponding formulae for a massless quark 
in the final state are presented in section |3|. In section [| QCD corrections to 
the angular distributions in charm and bottom decays are presented. For the 
sake of completeness we summarize our calculations based on [|l^ 
in Appendix A. 



and 23 



2 Angular and energy distributions 

In this section we give the formulae [] for the distributions of leptons from 
the weak decays of a polarised heavy quark Q of mass mi and the weak 
isospin ^3 = ±1/2. The mass of the quark q originating from the decay is 
denoted by and e = m 2 /mi. These formulae are written in the Q rest 
frame. However, the two variables which we use, namely the scaled energy 
x = 2£°/mi = 2Q ■ E/mf, where I is the four-momentum of the charged 
lepton, and the scaled effective mass squared of the leptons y — (£ + v) 2 /m\ 
are Lorentz invariants. For the massless lepton case considered in the present 
article also the third variable cos#, can be related to the scalar product i ■ s, 
where s^ 1 = (0, s) is the spin four-vector of the decaying quark and 9 denotes 
the angle between s and the direction of the charged lepton. S = \s\ = 1 

1 Fortran77 version of the formulae given in this article is available upon request from 
jezabek@chopin.ifj.edu.pl 
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corresponds to fully polarised, S = to unpolarised decaying quarks. 

The distribution of the neutrino for I 3 = ±1/2 is given by the formulae 
describing the distribution of the charged lepton from the decay of Q with 
the weak isospin I 3 = =Fl/2. 



2.1 Triple differential distributions 



The QCD corrected triple differential distribution for the semileptonic decay 
of the polarised quark with J 3 = ±1/2 can be written in the following way 



dr= 



G 2 m\ \V C km\> 



dx dy d cos 9 



32vr 3 (1 
2a x 



3n 



y/y) 2 + 7 2 \ 

Ff(x, y) + S cos 9 Jf(x, y) 



Fo(x,y) + S cos9J^(x,y) 



(1) 



Vckm denotes the element of the Cabbibo-Kobayashi-Maskawa quark mix- 
ing matrix corresponding to the decay channel under consideration. The 
effects of W propagator are included as a factor 1/[(1 — y/y) 2 + J 2 ], where 
y = m^/m 2 and 7 = T w /m w ; c.f. |T!J. These effects are important for the 
top quark fl9 , For charm and bottom quarks mi is much smaller than 
the mass of W boson and the four-fermion limit can be employed in which 
the factor mentioned above is replaced by 1. 

The Dalitz variables x and y fulfil kinematic constraints: 



0< x < x m = 1 — e 

0< y <y m = x(x m — x)/(l-x) 



(2) 



or 



0<y<(l-e)- 

W- < X < Wa 



(3) 



The functions w± and other useful kinematic functions are defined as follows: 

(l-y + e 2 )/2 



Po 

P3 
P± 



Po ± P3 = 1 



w-. 



Y p = I In (p+/p_) = In (p+/e) 
Y w = ±ln(w+/w_) = ln(w+/y/y) 
z m = (l-x)(l-y/x) 



(4) 
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The functions F^(x,y) and Jq (x, y) corresponding to Born approximation 
read: 



o k x > y) 



' o 



F n {x,y) 



J o (x,y) 



(x - y)(x m -x + y) 

(x - y)(x m - x + y - 2y/x) 



(5) 
(6) 
(7) 
(8) 



The first order QCD corrections and the corresponding functions F 1 (x,y) 
have been calculated in [TP| whereas Jf(x,y) has been given in [53]. In the 
course of the present work we rederived these old results and from them we 
obtained the much simpler expressions given in this article. The result for 
3i(x,y) is new. The formulae for Ff(x,y) and Ht(x, y) readQ: 



Ff(x,y) 



F±(x,y)% + Y / A^ n + Ai 

71=1 



(9) 



where 
$ 



$4 
$5 



5 

X] 

n=l 



2po 



Li 2 fi-lzfUii, fl-i^U L i 2 fl 



-Li, 1 



P+ 
1 - y/g' 
P- , 



P4 



(10) 



1 — x 



P- 



Li 2 - Li 2 (tu+) + 4Yplne 



+4 (l - ^Y p j In (z m - e 2 ) - 4 In z m 

Li 2 («?_) + Li 2 (w + ) - Li 2 (x) - Li 2 (y/x) 

Y P /p 3 
|lne 

|ln(l-x) 
|ln(l - y/ar) 



(11) 



2 Throughout this article polylogarithms are defined as real functions. In particular: 

U 2 {x) = - / dyln|l - y\/y 
Jo 

u 3{x) = / dyLi 2 (y)/y 
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Af = x + e 2 x 

At = -pUS + 2x + y + e 2 )-2e 2 y 

At = -(l-y)(3 + 2x + y) + 2e 2 (l + hx) + e 4 

A\ = 5 - 5x - e 2 (4 + 5x) - e 4 

At = -2xy + 9x- Ax 2 - 2y - y 2 - 7e 2 x 

A+ = {-Axy + 4y-y 2 + y 2 /x)/2 + e 2 [x + y-x/{l-y/x)]/2 (12) 
Bt = —x + e 2 x 

B+ = p\ [ (5 - 2x - y - 2y/x - 2/x) + e 2 (-l + 2/x) ] + e 2 (-xy - x 

+2y — y/x — y 2 /x) + e 4 (x + y/x) 
Bt = (l-y){5-2x-y-2y/x-2/x) + 2e 2 {-3 + 7x + 2/x) 
+t\l-2/x) 

Bt = -3 + x + 2/x + e 2 (A-9x-A/x) + e 4 (-l + 2/x) 

Bt = -2xy + 3x- Ax 2 + 6y - y 2 - 2y 2 /x - 7e 2 x 

Bt = {2-2x 2 + 2y -3y 2 -2y/x + 3y 2 /x)/2 + e 2 [-4 + x + y + 2y/x 

- x /(l- y / x )]/2 + e i (13) 

A\ = —2xy + x + y + y 2 + e 2 (x — y) 

A 2 = pU-5 + 2x-3y + e 2 )-2e 2 y 

Al = (l-y)(-5 + 2x-3y) + 6e 2 (l + x-2y)-e 4 

A~ = 5 + 4xy-5x + 3y-5y 2 -2y 2 /x + e 2 (-4-5x + lly)-e 4 

A~ = Qxy + 9x -Ax 2 -lly -2y 2 + 2y 2 /x + 7e 2 (-x + y) 

A 6 = (3xy + 2y-3y 2 -2y 2 /x)/2 + e 2 [3y-y/(l-x)]/2 (14) 

B^ = —2xy — x — 5y + y 2 — 2y 2 /x + e 2 (x — y) 

B 2 = pl[(3 + 10x + y + 10y/x-2/x) + e 2 ( y l + 2/x)]+e 2 ( y xy + x-2y 

+y/x + y 2 /x) - e 4 (x + y/x) 
B7 = (l-y)(3 + l0x + y + l0y/x-2/x) + 2e 2 (-l + 5x-4y-2y/x 

+2 fx) - e 4 (l + 2/x) 
Bl = -3 + 12xy + x -7y -y 2 -12y/x + 8y 2 /x + 2/x + e 2 (A-9x 

+7y - A/x) + e 4 (-l + 2/x) 
B~ = Qxy -9x-Ax 2 -y -2y 2 + 10y 2 /x + 7e 2 (-x + y) 
B 6 = (2-5xy-2x 2 + 2y + 7y 2 -2y/x-2y 2 /x)/2 + e 2 [-4-9y + 2y/x 

+y/(l-x)]/2 + e 4 (15) 
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2.2 Double differential x 



— 6 distributions 



For the top quark in the narrow W width limit 7 — > the mass squared 
of the leptons is fixed at value y = (m w /m t ) 2 . Thus, in this limit the triple 
differential distributions of subsection [O] reduce to the double differential x — 
9 distributions. In the four-fermion limit (m w — > 00) the double differential 
x — 9 distributions can be calculated from the following formulae: 



dr 



± 



dx d cos 9 



dy 



dr± 



dx dy d cos 9 32n 3 
2a 



G l m i , 12 / 

VcKM\ J( 



X 



-S cos 9 j (x) 



3n 



ft(x) + Scos9jf(x) 



where 



and 







xf/{l 


3o( x ) 






fo( x ) 




-x) 2 [3 


3o( x ) 


gX {x m 


-x) 2 [l 




ft 


(x) = 




3i 


(x) = 



2x + e 2 + 2e 2 /{l-x)}/{l-xf 



dy Ff(x,y) 
dy Jf(x,y) 



(16) 



(17) 
(18) 
(19) 
(20) 



(21) 
(22) 



The integrals in eqs.(|2~T|)-(|2"2"D can be calculated numerically. In the massless 



limit e — > analytic formulae for fi(x) and fi(x) have been derived. These 



formulae are listed in subsection |3.2| . The results for charm and bottom 
quarks have been already published in |15| |. For the polarisation independent 
parts of the distributions they are in perfect agreement with the results of 



rprgfl . Thus, they are in conflict with the results of |17| and (T^|; see section 5 
in [Tj| where a detailed comparison is given and cross checks descibed. 
A convenient way to present the results is to express the double differen- 
tial distributions as products of the energy distributions and the asymmetry 
functions, c.f. eqs.(5)-(7) in HI 



dr= 



dx d cos 9 



dT± 

dx 



[1 ± a±(x)cos#] /2 



(23) 



where 



a±(x) = ± S 



3oi x ) 



2 -Bjf( x )} I \fti x )- 2 -tft{ x 



(24) 



The functions a±{x) for charm and bottom quarks have been given in fT5|| . 
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2.3 Double differential y — 6 distributions 



We follow the normalisation convention of |2S[] and write the double differ- 
ential y — 9 distributions in the following way: 



dP 



dy d cos 9 



dx 



dr 



± 



\V ( 



CKM\ 



where 



and 



><2 [Mv) + S cos 9 Jt{y) 



dx dy d cos 9 
2a s 



3tt 



G 2 F rrv[ 
192tt 3 (1 -|//j/) 2 + 7 2 



F^ + Sco^Jtiy) 



MV) = 4p 3 [(l-e 2 ) 2 + y(l + e 2 

J+(y) = My) 

J -(y) = T {y)+4&y{yY w -pz) 



2y 2 



My) 
Jtiv) 



i A^T?± 



dx Ff (x, y) 



12/ dxJf(x,y) 



(25) 



(26) 

(27) 
(28) 



(29) 
(30) 



Eq.(|29|) provides a very non-trivial cross check of this calculation. It has 
been verified both analytically and numerically that for Ff(x,y) given by 
eq.(H|) the integrals in (^9j) are equal. Moreover, Ti{y) can be calculated in 
an easier way, see formula (2.6) in [p3| , if the phase space integrals over the 
momenta of leptons are performed first. This cross check is also fulfilled. It 
has been shown, see sect. 5 of [TBI, that when J-'iiy) is analytically continued 



it agrees with the analytic results of p(| for QCD corrections to the width of 
W boson. All these tests indicate that our formulae for Ff (x, y) are correct. 
Since all the calculations for both F 1 (x, y) and jf (x, y) were performed in 
parallel using the same algebraic computer programs we believe that the 
same is true also for jf (x,y). 



Analytic expression for Fiiy) has been first given in |P5| . For the sake of 
completeness we rewrite this expression using the notation of the present 
article and some identities for dilogarithms. Our formulae for Fi{y) and 
Jx{y) read: 



My) 
Ji + (y) 

Jr(y) 



My)^o + A^i + A 2 Y W + A 3 Y P + Ap 3 hie + A 5 p 3 (31) 
J+(y) + + B+Y w + BtY p + Bfp 3 \ne + B+Ps 

+ B+Y w Y p + B+Y w \ne (32) 
J ~(y) + B^x + B^Y W + B^Y P + B^p 3 \ne + B^p 3 

+ B^Y W Y P + BjY w lne + 4% 2 (B 8 + B 9 p /p 3 ) (33) 
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where 

= 2 [4Li 2 (2p 3 /p+) - AY p \n(2p 3 /p + ) - lnp_ In w+ + lnp+ In w_ ] p /ps 
+81n(2p 3 ) -21ny (34) 

= 4Li 2 (w_) - 4Li 2 (w+) (35) 

and 

■Ai = Fo(y)Po/P3 

A 2 = -8(1 -e 2 )[l + y- Ay 2 -e 2 (2-y) + e\ 

A 3 = -2[3 + 6y-21y 2 + 12y 3 -e 2 (l + 12y + 5y 2 ) + e 4 (ll + 2y)-e 6 

At = -12[l + 3y-4y 2 -e 2 (4-y) + 3e 4 " 

A 5 = -2 [5 + %-6t/ 2 -e 2 (22-9y) + 5e 4 l (36) 



Bf = 
Bt = 
Bt = 



5 - 9y 2 + Ay 3 - e 2 (8 - 8y + 6y 2 ) + e 4 + 2e 6 

1 + Ay - y 2 - e 2 (3 + 3y - Ay 2 ) + e 4 (3 - y) - e 6 

3(1 - y) 2 (l + Ay) + e 2 (ll - 5y 2 ) - e 4 (13 - 2y) - e 6 

(l-y)(l + Ay)-e 2 (A-y) + 3e 4 



-2 
-12 



15 - y + 2y 2 - e 2 (12 + 7y) - 3e 4 



Bt = -2A(l + y-e 2 )(2p 3 + ye 2 /p 3 ) 
Bt = -24(l- 2 / 2 -2e 2 + e 4 ) 



(37) 



Bl = 

B 2 = 

Bt = 

Bt = 

Bt = 

Bt = 

Bt = 
Bl 



B a 



= A 



5-A2y + A5y 2 + Ay 3 - 2e 2 (4 + 8y + 3y 2 ) + e 4 + 2e 6 
8 [-1 + 8y + 10y 2 + e 2 (3 - 9y - Ay 2 ) - e 4 (3 - y) + e 6 

2 [-3(1 - y)(3 - 23y - 4y 2 ) + e 2 (l + 2Ay + 5y 2 ) + e 4 (7 - 2y) + e 6 

-12 [3-23y-4y 2 -e 2 (6-y) + 3e 4 " 

-2 [-15 + 37y-2y 2 + e 2 (12 + 7y) + 3e 4 " 

-24 [2p 3 (l - 5y - e 2 ) - e 2 (l + y - e 2 )y/p 3 

-24 [(l- y )(l-52/)-2e 2 (l-y) + e 4 

Li 2 (^) - Li 2 (p_/ P+ ) ] + Li 3 ( W+ ) - Li 3 (w_) + 8Y p \n(w + ) 



-Y w [Li 2 (w-) + Li 2 {w + ) + 81n(2p 3 ) - 21ny - 41n(w+) ] 



Li 3 (l) + Li 3 



W—P— 

w +p+ 



- Li 3 (p_/p+) - Li 3 (w_/w + ) 



+2Y W [Li 2 (w + ] - Li 2 (w_) - 4Li 2 (2p 3 /p + ) + AY p ln(2p 3 /p + ) 

+2Y w \np + -AY p \nw + ] + AY p [Li 2 (^) - Li 2 (p_/p + ) ] (38) 
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3 Massless case 



3.1 Triple differential distributions 

For m 2 = the functions in eq. (|TJ) simplify considerably. The limit e 
can be easily obtained for Ff(x,y) and jf(x,y). In this limit 



and 



o K x >y) 
o fay) 
o fay) 



J (*,y) 



x(l — x) 

K&y) 

(x-y)(l -x + y) 
(x-y)(l-x + y-2y/x) 



(39) 
(40) 
(41) 
(42) 



F+ $ + x$i - (3 + 2x + y)$ 2 o3 + 5(1 - x)$ 4 + {-2xy 
+9x - Ax 2 - 2y - y 2 )$ 5 + y(A - Ax - y + y/x) /2 (43) 



Fi(x,y) = Fg % + (-2xy + x + y + y 2 )^ + (-5 + 2x - 3y)$ 



2o3 



+ (5 + Axy -5x + 3y- 5y 2 - 2y 2 /x)$ A + (<oxy + 9x - Ax 2 
-lly - 2y 2 + 2y 2 /x)$ 5 + y(2 + 3x - 3y - 2y/x)/2 (44) 

Jf{x,y) = J+<$> -x$ 1 + {5-2x-y-2y/x-2/x)<$> 2<>3 + (-3 + x 
+2/a;)$ 4 + {-2xy + 3x - 4a; 2 + 6y - y 2 - 2y 2 /x)$ 5 
+ (2 - 2x 2 + 2y - 3y 2 - 2y/x + 3y 2 /x)/2 (45) 



Ji {x,y) 



where 



J $ + (-2xy -x-hy + y 2 - 2y 2 /x)<& 1 + (3 + lOx + y 
+10y/x - 2/x)$2o3 + (-3 + I2xy + x - 7y - y 2 - 12y/x 
+8y 2 /x + 2/x)$ 4 + (6xy - 9x - Ax 2 - y - 2y 2 + 10y 2 /x)$ 5 



+ (2 - bxy - 2x 2 + 2y + 7y 2 - 2y/x - 2y 2 /x)/2 



$o = £ + 2Li 2 (x) + 2Li 2 (y/x) + ln 2 (i T $ 

$i = 2d + u 2 {y) - U 2 {x) - U 2 {y/x) 

$2.3 = \{l-y)Hl-y) 

$4 =| ln(l - x) 

$ 5 = \\n{l-y/x) 



(46) 



(47) 



Eqs.(lD-( 

in |i and eq.(3.3) in |3| . 



agree with those given in the literature; c.f. eqs.(3.9) and (4.9) 
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3.2 Massless limit for x — 6 distributions 



For e = explicit analytic formulae can be derived for the double differential 
x — 9 distributions. The functions fo{x), (x), fi{x) and jt(x) which 
appear in eq.(|16D are given by the following expressions, cf. eqs.(3.10) and 
(4.10) in |L9|: 



fo( x ) 
3o( x ) 
3o( x ) 



4> (x) 



x 2 {l — x) 
x 2 (3-2a;)/6 

fo + ( x ) 

x 2 (l -2x)/Q 
/ o + (x)0o(x) + (1 
+±x(10 - 19x) 



X 



±(5 + 8x + 8x 2 )ln(l -x) 



36x + 42x 2 - 16x 3 ) ln(l - x) 



+ ^x(82 - 153x + 86a; 2 



72 

3o( x 
i 



o(x) + (l-x) |(-3 + 12x + 8x 2 + 4/x)ln(l 



x 



+ 



12 ' 



- 2x - 15x\ 

x) + 4(11 -36x + 14x 2 



Jo K x m 
+ ^(-8 + 18x - 103x 2 + 78x 3 ) 



16x 3 - A/x) ln(l - x) 



2Li 2 (x) + 2vr 2 /3 + ln 2 (l 



x 



(48) 
(49) 
(50) 
(51) 

(52) 
(53) 

I 

(54) 

) 

(55) 
(56) 



3.3 Massless limit for y — distributions 

For e — > the functions in eq. (p5]) read, cf. eqs.(3.1)-(3.2) in [^5 



My) 

J + (y) 

J (y) 



2(i-z/) 2 (i 

Mv) 

2(i-y)(i 



2y) 



lly-2y 2 )-2Ay 2 lny 



(57) 
(58) 
(59) 



2^2 



7T 2 + 4Li 2 (y) + 21nyln(l - y) - (1 - y)(5 + 9y - 6y 



+ Ay(l -y- 2y l ) In y + 2(1 - y)\h + 4y) ln(l - y) 
J+(y) = - 2 ^ 2 (l-y)(l + y + Ay 2 ) + (l-y)(lh-y + 2y 2 ) 



(60) 



+ 4(1 



y)(5 



5y - 4y 2 ) Li 2 (y) + 8(1 -y)(2 + 2y- y 2 ) ln(l - y) hiy 



+ 8y(2 + y - y 2 ) Iny + 2(1 - y) 2 (5 + Ay) ln(l - y) 
Jr{y) = -l7r 2 (l + 6y + 9y 2 -Ay 3 ) + (l-y)(15-37y + 2y 2 



(61) 



5 - A2y + A5y 2 + Ay 3 + 18y 2 In y Li 2 (y) 



+ 2A0y 2 [Lj a (l)-Li 3 (2/)] + 4 
+ 8(1 - y)(2 - 13y - y 2 ) ln(l - y) \ny + 4y [4 - (19 + u 2 )y - 2y 2 
-2(1 -y)(l + 19y + Ay 2 ) \n(l-y) 



hi y 

(62) 
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4 Angular distributions 



Angular distributions are related to the double differential y — 9 distributions 



discussed in subsection 2.3 



dr± ^ a dT± («*s 

d y A„.A—a ( 63 ) 



d cos 9 Jo dyd cos 9 

with dr ± / dyd cos 9 given in eq . (p5|) . 

For the top quark and the narrow W width the double differential distribu- 
tions are proportional to Dirac delta function S (y — m^/mf) and the integral 
in ( |63| ) is trivial. For the charm and for the bottom quarks the four-fermion 
approximation can be used. 



Figure 1: a) QCD reduction of the total decay rate: the ratios TZ(e)/TZo(e) 
as functions of e for a s = 0.1, 0.2, 0.3 and 0.4 - solid, dashed, dotted and 
dash-dotted lines; b) QCD correction to the angular asymmetry ot~(e) for 
a s = 0, 0.2 and 0.4 - solid, dashed and dotted lines, respectively. 



In order to illustrate the size of the radiative corrections discussed in 
the present article we employ the four-fermion limit and write the resulting 
angular distributions in the following way: 

dr± 



dcos^ 



r ft(e) 1 ± cr^e) cos# 



(64) 



where T = G 2 p m\\ Vckm\ 2 /^84ir 3 . In Born approximation for completely 
polarised quarks (S = 1) 

TZo(e) = 1 - 8e 2 + 8e 6 - e 8 - 24e 4 lne (65) 

1 (66) 



a [e) 



1 - 12e 2 - 36e 4 + 44e 6 + 3e 8 - 24e 4 (3 + 2e 2 ) lne /Ko(e) (67) 
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QCD corrections decrease the total decay rate. In Fig. la the ratios 
lZ(e)/1Z (e) are shown as functions of e for a s = 0.1, 0.2, 0.3 and 0.4 as 
the solid, dashed, dotted and dash-dotted lines, respectively. The correc- 
tions to eq.flSop are of the order of one percent or smaller. In Fig. lb the 
function a~(e) is shown for a s = 0, 0.2 and 0.4 as the solid, dashed and 
dotted lines, respectively. 

A Derivations and formulae 

In this Appendix we sketch our derivation of eq. ([!]). The method has been 
described in our earlier papers jl9| and | 23l . Many formulae presented here 



have been already given in these articles. We rewrite them using the notation 
adopted in the present article. Some formulae are simplified and misprints 
corrected. In the following Q, q, £, v and G denote the scaled four-momenta of 
the decaying quark and the decay products: quark, charged lepton, neutrino 
and gluon. The common scaling factor is l/m-i, and Q 2 = 1, q 2 = e 2 , 
£ 2 = v 2 = G 2 = 0. Thus, the gluon is massless. However, at intermediate 
steps we introduce the scaled gluon mass Xg in order to regularize those 
expressions which are infrared divergent. We introduce also P = q + G, 
z = P 2 , s - the spin four- vector of the decaying quark and R = Q — s. Some 
formulae are given only for the weak isospin J3 = 1/2 of the decaying quark. 
This implies that the analogous formulae for J 3 = —1/2 can be obtained by 
replacing £ with v and vice versa. 

A.l Differential decay rate 

The QCD corrected differential rate is given by the following formula 

dr± = dr± + dr^ 3 + dr± 4 (68) 

where 

dr± = G 2 m\ \Vckm? M± 3 d1l 3 (Q; q, £, z/)/vr 5 (69) 
is Born approximation, 

dr± 3 = | a s G 2 m\ \V CKM \ 2 M% dTZ 3 (Q; q, £, z/)/tt 6 (70) 

describes the virtual gluon contribution and 

dTf A = | a s G 2 F m\ \V C km\ 2 Mf A dK 4 (Q; q, G, £, z/)/tt 7 (71) 

comes from real gluon emission. Lorentz invariant n-body phase space is 
defined as 

dK n (P; Pl ,p 2 , ..., Pn ) = 5^(P - Er) J] P 1 (72) 
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A. 2 Three-body contributions 

In Born approximation the rates for the decays into three-body final states 
are proportional to the expressions 



Mt, 3 

Mn 3 



q-vR-lj [{l-y/y) 2 + 1 2 
q-lR-v I \(l-y/y) 2 + 1 2 



(73) 



The formulae for F^(x,y) and Jq (x, y) can be easily derived from (|69|) and 
( |73|) when Dalitz parametrization of the three-body phase space 



d!Z 3 (Q] q, £, u) = -^dxdydad(cos9)dj3 



(74) 



is employed. Integration over Euler angle j3 is trivial, and non-trivial integrals 
over a read, see Appendix B in [f23|, 



(75) 



Jdas-q = 27cScos6[{l + y -q 2 )/2-y/x] 
J da s ■ v = 2irS cos6[y /x — (1 + y — q 2 — x)/2\ 



As a final step the scalar products in A^o3 are expressed in terms of the 
variables x, y and cos 9 



Q-s = 

Q.q = (l + g 2_ y)/2 
Q-is = (1-q 2 -x + y)/2 
Q-£ = x/2 



£. s = -xS cos 6/2 
£.q = (x-y)/2 
£-u = y/2 
v ■ q = (1 — q 2 — x)/2 



(76) 



and eqs.@-(||) follow immediately. 

The same steps are performed for the three-body radiative corrections which 
are given by the formula (|70|) , where 



M+ 3 = - [(1 -y/y) 2 + i 2 \ [q-v R-£H + e 2 Q ■ u R-£H + 
+ q-u(q-£ + R-lQ-q - Q-£R-q)H_ 
+ \e 2 (vl + Q-vR-t - R-vQ-i) (H + + H_)} (77) 



H 4 



Li 2 1- 



4(1- poYp/ps) In X G + (2p /p 3 ) 
-Li 2 (l - - Y P (Y P + 1) + 2(ln e + Y p ) (Y w + Y p ) 
+ [2p 3 Y p + (l-e 2 -2y) \ne}/y + 4 



l±(l-e 2 )/y\Y p /p 3 ±±\ne 



(78) 
(79) 
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A. 3 Four-body contributions 

Neglecting those terms whose contributions to the decay rates are 0(Xq In Xq) 
or smaller (i.e. vanishing in the limit A^ — > 0) one can cast the contribution 
of real radiation into the following expression 



Mt A - 



{ l_y/y)2 +1 2 



Bf Bt Bf 

1 ' ' (80) 



(Q ■ G) 2 Q-GPG (P-G)< 



where 



Bf = q-v[R-£(Q-G-l)+G-£R-Q-Q-£R-G + G-£R-G] 
Bf = q-u[G-£R-q-q-£R-G + R-£(q-G-Q-G-2q-Q)} 

+ R-£(Q-vq-G - G-vq-Q) 
Bf = R-£{G-v q-G - q 2 P -v) (81) 

The four-body phase space is decomposed as follows 

dft 4 (Q; q, G, t, v) = dzdTZ 3 (Q; P, £, v)Ml 2 {P- q, G) (82) 

After q is substituted by P — G integration of -M^ 4 over &R,2{P] q, G) is 
performed. Due to Lorentz invariance all the integrals which appear in this 
calculation can be reduced to the scalar integrals 

I n = J dn 2 (P;q,G)(Q-G) n (83) 



Explicit formulae are listed in Appendix C of [g3| for the reduction of ten- 
sor integrals / dK 2 (P; q, G)(Q ■ G) n G a , J dK 2 {P] q, G)(Q ■ G) n G a G> 3 and 
Jdn 2 (P; q, G)(Q ■ G) n G a G p G~< to the scalar ones. 

In the next step Dalitz parametrization is employed for the three-body phase 
space dlZ-iiQ; P, £, v) in the same way as described in the previous subsec- 



tion, but now for q replaced with P in eqs.( [74T )-(|76|). In particular q = e is 
replaced with P 2 = z. 

After integrations over Euler angles a and j3 the contribution of real radiation 
is splitted into two pieces: the infrared divergent part of the form 



Const [F (x, y) + S cos 9 J (x, y)] I div 

where 

Idiv = 1-2 ~ (1 - V + e 2 ) I-i/(P ■ G) + e 2 J /(P • Gf (84) 
and the rest which is infrared finite. 

For Ac ^ 1 the infrared divergent part can be integrated over z and the 
result is proportional to 
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dzl div = 4 1 - ™Y P ) In _ 2 In (z m /e 2 ) - 2 



7TJ(e+A ) V P3 J \ e ^G 

2p 



+ 



P3 



u,(i_iz£ + Li 2 (i-i^ -Li 2 fi-lz£ 
V P+ / V p+ / VP-, 

-n 2 ( 1 -i^)-L,( 1 -£-)- W + 1 , 



(85) 



When this result is added to the contribution of virtual corrections the in- 
frared divergent terms ~ fnAc cancel out. The sum is simplified using the 
following identity 

Li 2 (l-^)-Li 2 (l-^)-Li 2 (l-^) = 

Li 2 (tw_) - Li 2 (w+) - 2Y W \np + (86) 

For the infrared finite part the limit — > is performed. The formu- 
lae for the scalar integrals I n simplify considerably in this limit and read 
1-2 =tv/(P-G) , ' I^ = 7ry p (z)/V^A , 

T APG) n+l \ (Q-P + v^A) n+1 -{Q-P- v^A)" +1 ] . >m 
In= 2(n + l)(P^v^=A (n - 0) (8?) 



wheref] 



A = Q 2 P 2 - (Q ■ P) 2 = -A(z)/4 
\{z) = z 2 -2(l + y)z+(l-y) 2 



y P (z) = |ln 



l-y+z+W\(z) 



P-G = (z-e 2 )/2 (88) 

Integration of the infrared finite part over z is tedious. It would be difficult 
to accomplish this task without FORM |27|. All the integrals which appear 
can be reduced to integrals / dz z n / X m (z) and / dz z n y p (z) / X rn+1 ^ 2 (z). Thus, 



they can be derived from the recursion relations given in Appendix B of [IE . 



A. 4 Formulae for Ff(x,y) and jf(x,y) 

After some algebra mentioned in preceding subsections one derives the fol- 
lowing formulae for F 1 (x, y): 

Ff{x, y) = Hf(x, y) + W±(ar, y, e 2 ) - Hf(x, y, z m ) (89) 



3 In the present calculation only n < 1 is needed. Thus, instead of eq.(^7|) one can use: 
I„ = n(P ■ G) n+1 (Q ■ p)«/(p 2 )«+! for (n = 0, 1) . 
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where the first term in r.h.s. is the sum of the three-body virtual corrections 
and the infrared divergent four-body ones. The two other terms originate 
from the infrared finite four-body piece. For jf(x,y) one has: 



Jj (x, y) = K x (x, y) + K, 2 (x, y,e)-)C 2 (x, y, z m ) 



(90) 



The results of |19| and [g3[] are presented in this way; see eqs.(3.2)-(3.6) and 

(4.5)-(4.8) in and (3.1),(A.1)-(A.15) in J£§. 

Both H 2 (x, y, z) and K 2 (x, y, z) are given by lenghty expressions: 

nt(x,y,z) = (af + zai)y p (z)/\V 2 (z) + (af + zaf) y p (z) / \f\(zj 

+ afy p {z)^X(7)+ (4 + zaf)/\{z) 
+ a±z + a±z- 1 + a± In z + a± [Li 2 (W + (z)) + Li 2 (W_(z))] (91) 



1Cf(x, y, z) = (bf + zbf) y p (z)/\ 5 / 2 (z) + (bf + zbf) y p (z)/\ 3 / 2 (; 



+ (bf + zbi)y p (z)/^x(z)+ b±y p (z)y/xw 

+(bf + zbf)/X 2 (z) + (bf + zbf 1 )/X(z) 
+ bf 2 z + bf 3 z- 1 + bf 4 lnz + bf 5 [U 2 (W+(z)) + Li 2 (W_(z))] 

where af and a, are rational functions of x, y and e 2 , and 



(92) 



WUz) 



l + y-z±JX(z) 



/2 



(93) 



These complicated expressions simplify dramatically for z = e 2 and z = z m . 
At these special points 



^(e 



p 

W±(e 2 ) = w± 
X(z m ) = (x- y/x) 2 

y p (z m )/ \l X(z m ) = 2^y/.r) 

Li 2 (W+(z m )) + Li 2 (W_(z m )) = Li 2 (x) + U 2 (y/x) 



(94) 



and after some algebra eqs.(||) and (|10f) are derived. 
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